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Abstract 

A new approach leading to the formulation of the Hamilton-Jacobi 
equation for field theories is investigated within the framework of jet- 
bundles and multi-symplectic manifolds. An algorithm associating 
classes of solutions to given sets of boundary conditions of the field 
equations is provided. The paper also puts into evidence the intrinsic 
limits of the Hamilton-Jacobi method as an algorithm to determine 
families of solutions of the field equations, showing how the choice of 
the boundary data is often limited by compatibility conditions. 
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I Introduction 

The Hamilton-Jacobi equation in classical point particle mechanics is a pow- 
erful tool, used to solve the evolution equations, by means of the concept of 
complete integral. Moreover, its theoretical importance is placed in its being 
the cornerstone of the transition from classical to non-relativistic quantum 
mechanics in Schroedinger' s formulation. 

An analogue of Hamilton-Jacobi equation for field theory has been cre- 
ated at the beginning of the XIX century through different alternative ap- 
proaches In particular the formulation of De Bonder and the one of 
Caratheodory have been rather studied and further developed, even recently. 
Yet none of them has proved to be as powerful as the theory which is gener- 
ally used in mechanics. In fact, while the equation itself seems to be even less 



constraining than the mechanical one, being a single PDE in many unknown 
functions, an additional set of embeddability conditions is needed to relate its 
solutions to the field equations. This results into a set of quite restricting 
conditions and the possibility of finding solutions strongly depends on the 
nature of the particular problem under consideration. 

It is clear in the literature that the construction of a parametric set of 
solutions, resembling the classical complete integral, is quite often impossi- 
ble, due to the presence of strong integrability conditions. Therefore, many 
different authors pursued the aim of producing sets of particular solutions: 
some of them built single solutions starting from given solutions of the field 
equations ^ ^, others use the concept of Baecklund transformation under 
particular dimensional restrictions 

The present paper is aimed at building a particular class of solutions of 
the Hamilton-Jacobi equation, extending the previous approaches. What we 
show here is an algorithm associating particular solutions of the Hamilton- 
Jacobi equation to given sets of boundary conditions of the field equations. 
This approach differs from the previous ones in its being independent of any 
choice of the possible surfaces where the boundary data are assigned and of 
every dimensional consideration. Moreover there is no need of knowing any 
particular solution of the field equations in advance. 

Unluckily, the argument is not entirely straightforward, since the above 
mentioned embeddability conditions play a very strong role and eventually 
put bounds on the possible initial data which are able to generate solutions 
of the field equations by way of Hamilton-Jacobi. These topics are examined 
within the proposed class of solutions, which is by no means the most general 
one. It is opinion of the author that some generalizations of the algorithm 
could be possible and could even help weakening some of the restrictions; yet 
most of the limitations do not actually depend on the particular choices, but 
on the nature of the problem itself. 

In section 2, the Lagrangian and Hamiltonian aspects of field theories 
within the framework of jet-bundles are revisecP'^'^. These arguments are an 
actual survey on currently known arguments and are reported for notational 
purposes onlji^. 

Section 3 is devoted to the study of the Hamilton-Jacobi equation. After 
proposing an easy alternative deduction of the equation itself, completely 
equivalent to De Bonder's one, the concept of embeddability of solutions is 
explained and the conditions for this to happen are worked out. Then, the 
algorithm leading to the determination of a class of solutions is explained in 
detail. 

Finally, an easy example is provided in Section 4, with the aim of showing 
how the method works, avoiding useless computational complications. 
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II Preliminaries 



This section is devoted to recalling the geometrical environment where La- 
grangian and Hamiltonian dynamics for field theories is developed, within 
the framework provided by jet-bundles. 

(i) It is known that field theories are generally described from a geometric 
viewpoint taking a bundle tt : E —>■ M into account, endowed with a set of 
local coordinates x^, . . . , x", y^, . . . , y"^, the first set denoting the basis coor- 
dinates on M and the remaining ones the local coordinates on the fibers. A 
physical field is then represented as a section if : M ^ E, locally written in 
the form: 

f = ^\x^') i = l,...,r, /i = l,...,n. 

(ii) Lagrangian theories are described taking the first jet-bundle ji{E) as- 
sociated with 7C : E —>■ M into account and working on the geometrical 
framework provided by the double fibration ji{E) E M. For varia- 
tional purposes, the manifold M is supposed to be endowed with a metric 
= Qfiudx'^ ® dx'^. 

The bundle ji{E) is endowed with a set of local fibered coordinates x^, y\ y^; 
the last set of coordinates represents the equivalence classes of sections ip 
having a first-order contact. The whole set is subject to the following trans- 
formation laws: 

x^ = x^ixn , f = y\x^,y^) , yl = ( + ^y^^ (1) 



dy^ ^ ) dx 

Eqs. (IH) show that the bundle ji{E) — »• E has the nature of an affine bundle, 
modelled on T*{M)(^V{E), where V{E) denotes the vector bundle of vertical 
vector fields on E —>■ M. 

The bundle ji (E) is endowed with a set of canonical differential forms, hence- 
forth called contact 1-forms, that result to be invariant under a coordinate 
transformation as the one described by eq. ([T]). They can be locally described 
in the form: 

uj' = dy' - yl dx^" . (2) 
Every section ip : M ^ E can be raised to a section ji (y?) : M ji {E) as 



dxf^ 



Conversely, every section (p : M ^ is said to be admissible if there 

exists a section ip : M —>■ E such that = ji(v?)- It is easy to prove that a 
section is admissible if and only if 0*{uj'^) = 
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(iii) Lagrangian field theories are described assigning a function y*, y^) 
on ji{E), called the Lagrangian, accounting for the dynamical aspects. De- 
pending on the theory under consideration, the choice of the Lagrangian is 
generally made following geometrical guidelines, symmetry properties or a 
constitutive analysis of the interactions. In this work we will sistematically 
suppose this choice has been made, and we will be only interested in its 
mathematical consequences. 

Equally important is the presence of a canonical n-form 6^ on ji{E), called 
the Poincare-Cartan form, induced by the presence of the Lagrangian func- 
tion and locally represented as 



dL 

6l = Lds + 7-— cj* A ds^ 



(3) 



where ds := y/g dx^ A ... A dx'^ is the volume form on M and cis^ := J ds. 
The latter is invariant under arbitrary coordinate transformations ([T]). Ac- 
tually, its importance lies in the fact that the variational principle for the 
functional I[ip] = J ^*{6i) on ii{E) singles out all the admissible sections 
if : M ^ E satisfying the Euler-Lagrange equations 



d 
dx^^ 



dL 



dy\ 



/!/ J 



dL_ 







(4) 



They have the nature of a set of second order partial differential equations 
for the functions denoting the local representatives of the section ip. 

(iv) Hamiltonian field theory is developed within the framework provided by 
the bundle K^{E) of skew-simmetric n-forms on E (compare withl^. Let us 
take the bundle A"(£') of horizontal n-forms into account, whose elements 
are annihilated whenever one of its arguments is vertical. Then, let A2(-E') 
be the bundle of 1-contact forms on E, whose elements have the property 
of vanishing when at least two of their arguments are vertical vector fields. 
The elements of the two bundles can be locally represented by the following 
conditions: 

a G A"(i?) ^ a = po{a) ds , (5a) 

a e A^iE) ^ cr = p{a) ds + (a) dy' A ds^, . (5b) 

The bundles A'^{E) and A2 (-E) are respectively described by a system of local 
coordinates x'^, y\po and x^,y\p,p^, subject to the following transformation 
laws: 

d{x^, . . . , x") 



Po = Jpo 



where J := det 



(9(a;^, 







p = J [ p + p 



_p dy^ dx^ 



^ dx-^ dxP 



Pi 



J dy^ dx^ _p 
dy"- dxP^^ 



(6a) 



(6b) 
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The latter shows that A"(i?) E is a. vector sub-bundle of A2{E) E; the 
quotient bundle Il{E) := A2 (-E')/A"(£') will be henceforth called the phase 
space and endowed with a set of local coordinates x^,y\p^. The transfor- 
mation laws (lEb) make A'2(-E) ^{E) into an affine bundle, modelled on 

(v) The bundle Ag is endowed with a canonical Liouville n-form 0, locally 
expressed as 

e =pds+p'^dy' Ads^ , (7) 
whose differential is represented by 

Q = de = dp Ads + dpf^ A dy' A ds^ . (8) 

The latter is a multi-simplectic {n + l)-form on A2{E) endowing it with a 
non-singular multi-symplectic structure. 

(vi) The presence of the Liouville n-form allows to write the field equations 
(jlj) on A2{E). In fact, there exists a unique fibered differentiable application 
A : ji{E) A2 such that A* (9) = 9^. This requirement can be translated 
into local fibered coordinates, representing the map A in the form: 

The application A : ji{E) —>■ A2{E) represents a sub-manifold of A2{E), 
fibered on the phase space n(-E); its image will be denoted by S := X{ji{E)). 
In particular, whenever the regularity condition 



det 



dyidy"^ 



+ (10) 



is satisfied, the first equation ([9]) can be locally inverted, allowing to obtain 
the coordinates = yl^{x'^, y\Pi)- If the latter is substituted into the second 
equation (Q, the following cartesian representation for the manifold S is 
obtained: 

S: p = L{x^,y\y;{x\y\p>;))-y;{x'',y\p';)p>; = -H{x^,y\p^) . (11) 

The function if(x'^, ) will be called the Hamiltonian of the system: it 
actually is (the pull-back of) a function on Il{E). 

A comparison with eq. ffTTj) shows that whenever the regularity condition 
ffTUj) is fulfilled, the sub-manifold S has the nature of a local section of the 
affine bundle Al^^E) —> Il{E), locally represented as: 

: n(E) A^(E) ; h:p + Hix'',y\p':)=0 . (12) 
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The latter will be henceforth called the Hamiltonian section and its assign- 
ment is completely equivalent to the knowledge of the Lagrangian function 
on ji{E), because of eqs. ([9]) and (ITTll . 

(vii) Finally, the presence of the Hamiltonian section allows to pull-back the 
Liouville 1-form on Il{E) and to build a variational principle on the phase 
space. Consider a section ip : M —>■ Il{E), locally represented in the form 

^: f = y%xn , p': = P^{xn . (13) 
An easy calculation shows that the stationary sections for the functional 



(14) 



are the solutions of the following differential equations 

dy' _ dH dpf _ dH 



which are the equivalent of Hamilton equations in field theories. They have 
the nature of a system of first order PDEs in the unknowns y^ and p^ and 
are often harder to solve than their Lagrangian counterpart, because of the 
difficulties needed to uncouple them. 

Every solution tp of the sistem (fT^ is A-related to the solutions of eqs. 

This means that the determination of a solution of ffT^ is equivalent to 
determining a solution of eq. Every solution (p : M n(i?) of eqs. ffT^ 
will be called a critical field. The determination of an alternative method 
yielding solutions of ( |T^ is the principal aim of the present paper. 



Ill The Hamilton-Jacobi equation 

The usual deduction of the Hamilton-Jacobi equation in analytical mechanics 
is based on the use of canonical transformations on symplectic manifolds 
by means of the concept of generating function. Unluckily, no such theory 
is available in the multisymplectic case. Our deduction of the Hamilton- 
Jacobi equation will be therefore based on an alternative approach, relying 
on the possibility of determining a particular sub-manifold of the surface 
iS, containing the lift of a solution of the field equations with the given 
boundary conditions. This requirement singles out a first order PDE and a 
set of additional embeddability conditions. 
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Let us consider the geometric structures determined by the following con- 
struction. Let cr G A"^^(£') be a generic horizontal {n — l)-form on E, 
locally described in the form: 

a = S''{x^',y')dsa • (15) 

Its exterior differential da G A2 (-E) is locally described as: 

da = — dsH — dy Adsa ■ (16) 

ox°' oy^ 

It represents a section : E —* A2{E) (compare with (IHb)) locally written 
as: 

Such a section is the natural generalization of the concept of Lagrangian sub- 
manifold to the multi-symplectic case, since the pull-back of the Liouville 
n-form to S is exact: 

We are now going to determine under what circumstances the above section 

: E ^ A2{E) is contained in the sub-manifold S C A^iE). 
A comparison between eq. (fTTl) . defining the sub-manifold S, and eq. ([1] 
shows that the following condition must necessarily be satisfied: 

+ H{x'',y\—-] =0 . 



dx°' \ dy 

The latter is known in the literature as Hamilton- J acohi equation for field 
theory. This equation is deeply linked with the solutions of the field equa- 
tions; in fact, performing a derivative of both sides of (fTSj) with respect to 
we obtain that: 

d dS^ OH d^S" dH d fdS'' , „ ,„A dH 

- + ^+ ^ . ^ = ^ ^ ^r-{x^. yW) 



dx"" dy^ dy' dy'-dy^ dp"^ dx"" \ dy"- J dy 



This means that the quantities = fulfill the second set of eqs.f ll4p for 
every y^ = y^{x^). Substituting them into the first set of ([T3j) we obtain the 
condition: 



Q^fl Qpt^ \y ' ' Qyi 

Eq. fliyp is not identically satisfied by all the solutions of the Hamilton-Jacobi 
equation, but is needed in order to make the field equations valid. 
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Definition 1 A solution = y^i^x"), = Pi{x^) is said to be embeddable 
in a solution = S^{x^,y^) of (ITSl) if condition (fT^ holds. 

If /i = 1, i.e. in tlie case of point particle meclianics, eq. (fT^ describes a veloc- 
ity field on E for every solution of (ITSl) and its integrability is guaranteed by 
Cauchy theorem; therefore, it is possible to determine a family of solutions 
of the field equations embedded in every solution of the Hamilton-Jacobi 
equation, at least locally. In the general case, instead, different scenarios are 
possible. 

1. If the integrability conditions for eq. ( |T9l) are satisfied, then a family of 
solutions of the field equations can be embedded in a solution S'^{x^,y''). 
This condition is called in the literature strong embeddability condition and 
is generally satisfied in much limited situations only; however this condition 
is absolutely required to build a complete analogue with the point-particle 
mechanics, i.e. to determine a complete integral of the Hamilton-Jacobi the- 
ory. This aspects of the theory will not be analyzed here. 

2. In general, the integrability conditions are not satisfied and this gives rise 
to a restriction on the set of possible solutions. In particular, the generic so- 
lution of the Hamilton-Jacobi equation is not able to single out a parametric 
set of solutions matching all the possible boundary data. This limitation is 
intrinsic with the nature of the equation itself and an analysis of the possible 
solutions needs to be handled case by case. 

3. If a generic solution of the Hamilton-Jacobi equation is given, condition 
ffT^ can be checked on itP; this last possibility is called weak embeddability 
condition and ensures the validity of the field equations only on the ac- 
counted solution ?/* = y^{x^). The problem with this approach is that if the 
integrability conditions for eq. ( |T9l) are not satisfied around the given critical 
field, there is no way of using the obtained solution of the Hamilton-Jacobi 
equation to generate other critical fields. Therefore, this approach has no al- 
gorithmic application, since solving the Hamilton-Jacobi equation becomes 
useless if its only embedded solution is already known. 

The approach adopted here will follow an intermediate viewpoint: we will 
determine a parametric set of solutions of the Hamilton-Jacobi equation, 
associated to a set of given boundary conditions. This is made possible by 
the fact that eq. ( |T8|) has the nature of a single first-order partial differential 
equation in the n unknown functions S°'{x'^,y^) and therefore possesses a 
great amount of equivalent solutions. This arbitrariness will be eventually 
fixed by the weak embeddability condition. Since only the solution associated 
with the given boundary data is sought, no integrability condition for fll9p is 
required. 
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In order to associate a solution to a set of boundary conditions, we first 
remind that they essentially consist in the assignment of the fields and of 
their normal derivative on an {n — l)-dimensional submanifold i : N —>■ M, 
locally expressed as 

where . . . , (""'^ represent a system of local coordinates on A^. The whole 
set of boundary data is therefore represented as: 



(20) 



x'' = ip''{C^) A = l...n-1 

y' = i^'{C^) /x = l...n 

where n = n^-^ is a vector field defined on and transverse to it. 
Now, let us consider the particular solution of the Hamilton-Jacobi equation 
obtained through the following procedure. Let X = X^{x^)^^ e D^[M) 
be an arbitrary vector field having the property of being transverse to the 
surface of the boundary data (but not necessarily equal to n on N) and take 
the following class of candidate solutions of the Hamilton-Jacobi equation 
into account 

S^{x'',y^):=^{x'',y^)X'^{x'') , (21) 

where (p{x^, y^) is an arbitrary function on E. Eq. (I2T1) fuUfills the Hamilton- 
Jacobi equation if the function (p{x'^, y*) is such that 

Eq. (1221) leads to a great simplification of the problem; in fact, it has the 
nature of a single first order partial differential equation in the single unknown 
ip{x^, y*) and can be easily solved using the same techniques as the traditional 
Hamilton-Jacobi equation (e.g. the separation of variables), for every choice 
of the vector field X. 

We will now show how determining the solutions of eq. ( l22l) can be useful 
to evaluate the critical fields. For this purpose, we will develop a solution 
algorithm, following the steps below: 

• we think of the components of the vector field X^ as free parameters and 
determine a set of solutions = ip{x'^,yi,X^), depending on arbitrary inte- 
gration constants; 

• we search for (at least) a vector field X^ = X^{x'^) satisfying the weak 
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embeddability condition on the given solution; this procedure is partly de- 
pendent on the choice of the boundary conditions, since the vector field X 
needs to be transverse to the surface A^; 

• we calculate the integration constants, imposing that the solution matches 
the boundary conditions. 

The parametric role played by the quantities suggests that they should 
appear algehrically in eq. ( l22l) . Therefore we choose to impose the condition 

on the vector field X; it makes the weak embeddability conditions harder, 
but simplifies the analysis and the solution eq. ( !22l) . 

We will cope with the Hamilton-Jacobi equation using the method of char- 
acteristics for first-order partial differential equations, stating the argument 
from a geometric viewpointP. Let us take the manifold T{E) into account: 
the latter is endowed with a set of local coordinates x'^,y\ X^ and these are 
the actual variables appearing in eq. (1^ . In fact, the choice provided by 
eq. (l23l) avoids the necessity of using an additional set of variables Xj^. 
Let P — *• T{E) denote a bundle of scalar functions over T{E) and let u be 
the coordinate along the fiber. Every solution of the differential equation can 
be thought as a section ip : T{E) P. Once the choice of the base manifold 
T{E) has been made, the presence of eq. fl2^ also avoids the necessity of con- 
sidering the quantities among the unknowns of the differential equation: 
this would actually bring us back to the original problem of having multiple 
unknowns. 

Provided the above choices are made, the first jet-bundle ji{P) is then en- 
dowed with a set of local fibered coordinates x^,y\X'^,u,Ufj_,Ui,U^, with 
the formal identifications ~ ^'""^ — ^'^m — "^^^^ is a classical 

context where PDEs can be framed from a geometrical viewpoint, interpret- 
ing the differential equation as a sub-manifold of ji{P), locally written in 
cartesian form as 

y\ X^ n, n^, u,, U^) = u^X^ + E{x^ , y\ X'^) = . (24) 

We notice that eq. ([21]) independent of the variables f/^. The characteristic 
curves can be calculated solving the usual system of first order ordinary 
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differential equations on ji(-P): 

<h^_dF_ _ dy^_dF_ _ dX^ _ dF 

du^ ' dui ' d^ dU^ 

du dF dF dF dF 
— - — + ^i^ir- + u,— + U,,- 



d^ du dUfj, dui dU^ 

dM^__dF_ d£ _ d/a^ _ _dF^ _dF_ _ dU^ _ dF dF 

Iti ~ ~d^ ~ 'd^^^ ' lli ~ ~ 9^^' ' 'd^ ~ ~dX^ ~~d^ 

(25) 

In tlie present case, taking eqs. and (1^ into account, tlie above equa- 
tions can be evaluated as follows: 

dun dH du dH 

Tlie first equations (!2Ub ) allows to determine the vector fields X'^, which 
result to be constant along the characteristics, and show that the projection 
on M of the characteristic curves are represented by their integral lines. 

X^ = A^{C'') x^ = A^{C^)^ + x^,{C^) (27) 

The quantities Xq^i^"^) can be easily evaluated by imposing ^ = in ( 1271) and 
comparing with the boundary conditions (l20l) . Moreover, it is also possilble 
to locally determine the values of ^ and as functions of x^ whenever the 
condition 

/ A^(C^) \ 

det dAH.^dxl ^0 (28) 
V dC^^^ dC^ I 

is satisfied. Being the quantities C,^ independent coordinates on N , it is easy 
to prove that eq. ( l28i) is satisfied whenever the vector field X is transverse 
to the sub-manifold N itself. 

Substituting eq. (128|) into (l26b) we obtain a system of 2n ordinary differen- 
tial equations of the first order in the unknowns y*(^, C"^, X^), Mj(^, X'*), 
whose local solvability is guaranteed by Cauchy theorem. Notice that the 
variables C"^,X^ only play a parametric role, being the latter ones indepen- 
dent of ^. The solution is provided by the families of functions 

l/^=f(e,c^^^«„/30 
«.=^?.(e,c^x^«.,/50 
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written in terms on 2n integration constants ctj, We now possess all 
the elements that are needed to impose the embeddability conditions (1191) . 
Taking the regularity condition (!28l) into account, they are equivalent to 
imposing 

Qyi gr^^l Qj^ Qj^fi Qj^ 



which is identically satisfied because of (l26b) and 



, (30a) 



that needs to be imposed by hand on the solution (!^ . Notice that the latter 
is already dependent on the choice of the boundary data. 
As a matter of fact, eqs. (130b) are a set of r x (n — 1) equations and the exis- 
tence of a vector field X'^ = X^{(^) fulfilling them is not guaranteed. If this 
is not the solution may still be determined adding some of the integra- 

tion constants Oj, /3* among the unknowns. We remark that this drawback is 
intrinsic in the nature of the problem, since condition ffTIJl) does not depend 
on the particular algorithm we use, but on the possibility of embedding the 
critical fields into the solution of the Hamilton- Jacobi equation. On the other 
hand, the translation of eq. ( fT9l) into (l30b) depends on the particular choice 
represented by eq. (!22|) : the existence of other possible algorithms leading to 
a weaker embeddability condition cannot be excluded. 

In the follow up, we will suppose that a vector field satisfying (IHUb) exists, 
though not excluding that its existence imposes an eventual restriction on 
the integration constants. Then, the following proposition holds: 

Proposition 2 Let X^ = X^{x'') be a vector field satisfying eq. ( !23ll . For 

every choice of X^ satisfying eq. (Wb). and for every ai, (eventually re- 
stricting to the ones making eq. (2Fb) solvable), the functions 

y' = f{axn,c^{xn,x^{c^{xn),a,,(3') 

are a solution of the field equations. 

The remaining problem is that of matching the above solution with the given 
boundary conditions (1201) . If no restriction on the possible choices of the 
integration constants at, /5* is present, they can be determined for every choice 
of the boundary data on the given surface. Otherwise, eq. flHUb) singles out 
a compatibility condition on the boundary data for each given solution and 
in some cases the problem could even turn out to be impossible. Up to now, 
this situation has still to be handled case by case. 
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The solution of the Hamilton- Jacobi can be obtained integrating eqs. fl^Ub ): 
the first equation (pUb ) allows to determine the quantities = u^{^, 7^), 
associated with the given solution and the given boundary conditions, the 
quantities 7^ being arbitrary integration constants. 

Substituting into the second equation (l26b) we can obtain the solution of the 
Hamilton-Jacobi equation evaluating the integral: 

«(e,C) = J [u.X'^ + u^^X'^^ di . (31) 

Taking eqs. ( |20l) and ( |271) into account, we can write the solution as 

= u{i{x'')X\x^)) X^{C\x'')) . (32) 

We remark that eqs. (l26l) are not linear in the quantities X^; therefore no 
superposition is possible among different solutions obtained using different 
fields satisfying the embeddability conditions. 



This paper suggests a method to find a certain class of solutions of the 
Hamilton-Jacobi equation, extending the ones that are currently available 
in the literature; as far as the present status of the topic is concerned, this 
does not represent the most general approach possible. Further developments 
could be eventually achieved modifying the form of the given solution (I2TII 
in its dependence in the (p and X^. Unluckily, every tentative beyond the 
linear dependence in the complicates the relationship between the char- 
acteristics and the integral lines of the vector field X: the problem still needs 
to be analyzed in detail. 



IV An example 

We present here the simplest example possible, in order to show how the 
method works and what kind of solutions can obtained, avoiding unnecessary 
techical complications. For this reason, the example is limited to dim{M) = 2 
and the metric is taken to be flat, diagonal and positive definite. None of 
the above simplifications is needed to make the algorithm work. 
Consider a free scalar field y = ?/(x^,x^), whose Hamiltonian can be written 
as 

^f(x\x2,i/,p\p2) = i((5^,pV + /iV) , /^=1,2 . 
The Hamilton-Jacobi equation takes the form 
dS^" 1 dS^dS 



dxf^^ 2[^'"' dy dy 
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Taking eq. (!?T|) into account, the latter can be written as 



This last equation could be solved through the separation of variables. In- 
stead, we will follow the guidelines provided in the paper. 
Let us suppose the boundary data are given on the surface x^{z) = 0, x^{z) = 
z in the form: 

y{0,z) = f{z) ; ^{0,z)=g{z) . (33) 

Taking the boundary conditions as well as eqs. ([27]) into account we have that 
the projection on M of the characteristic curves can be locally expressed as 
functions of the components and of the vector field X as: 

x\z,0=A\z)^ ; x'iz,0=A'iz)^ + z . (34) 

This last equation can be used to locally invert the Jacobian matrix | ^g'^g'^-^ ^ | 
and to calculate the partial derivatives with respect to C,,z. This allows to 
impose condition ( l23l) on the field X as follows: 

j^i^^^^^M = A\z)^^^ A\z) = e''A\z) . (35) 

We can now write eqs. fl26b ) explicitly: 

They give rise to the following sets of solutions: 

j/(^, z) = a{z) cos(/i a;(z) ^) + h{z) sin(/i a{z) ^ (36a) 



Mi(^, z) = —7— {b{z) cos{jja{z) ^) — a{z) sm{jj,a{z) ^)) (36b) 
a{z) 



where 



The weak embeddability condition gives rise to the following system: 

aiz)^^^ = A\z)^^^^il + e'''') , (37a) 

-b{z) , (37b) 



da{z) 




dz 


dz 


da{z) 


= e'A\z)-^. 


dz 


a{z) 
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db{z) 



fx 



■a{z) 



(37c) 



dz ^ a{z) 

Eq. (l37k ) puts a relationship between a(z) and ^^(-z) and can be easily 
integrated as follows: 



a{z) 



,2c 



Substituting into (IHTb.c) we obtain the following set of differential equations 
for a{z) and b{z): 

da{z) jjie'^ db{z) 



dz 



fie 

VTTe- 



2c 



b{z) 



dz 



fi e 

TrTe 



2c 



a{z) 



(3^ 



These equations are independent of the choice of a{z). Therefore, eqs.(l38l) do 
not put any restriction on the quantities a{z) (that are related to the vector 
field X), but result into a set of compatibility conditions on the boundary 



data. In fact, solving eqs. 
only if 



we obtain that 0361) satisfies the field equations 



a{z) 



b{z) 



A cos 



B cos 



/ie 



2c 



,2c 



z ] + i? sin 



2c 



VTT 



,2c 



/ie 



2c 



VT+e 



2c 



A sin 



/ie 



2c 



VT+e 



2c 



where A and B are arbitrary constants. This means that, evaluating for 
^ = 0, the only possible boundary data are of the form 

f{z)=a{z) ; g{z) = fib{z) . 

Substituting into eqs. ( j36l) we obtain the most general solutions of the field 
equations with the given initial data, which are derivable from the Hamilton- 
Jacobi approach. Even in this simple case, they are by no way the most 
general ones. 

Integrating eqs. fpUb ) we obtain that u^{z,^) = 7^ = const.; then the solution 
of the Hamilton-Jacobi equation can be obtained evaluating the integral 

After the rather tedious process of integration, the solution can be obtained 
taking eqs. (IMl) into account and evaluating functions oi z,^ for any 

possible choice of the function a{z). We obtain that: 



u{z{x^,x^),i{x^,x^)) 



,2c 



S'^{x^^x^) = u{z{x^ , x^) , ^{x^ ^ x^)) 



VTT 



,2c 
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